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APPLICATION OF A NEW MODIFIED PARAMETRIC ALGEBRAIC MODEL FOR 
MAGNETORHEOLOGICAL DAMPER IN VEHICLE SEMI-ACTIVE SUSPENSION 

SYSTEM 



The purpose of this paper is to present a new modified parametric algebraic model for MR damper in 



controlling the semi-active suspension system. Except for its accuracy, modified algebraic model is additionally more 
preferable in terms of its low computational expenses compared to differential modified Bouc- Wen's model that is 
extremely computationally demanding. The control part consists of two nuzzled controllers, one the system controller 
that generates the required damping force, and other the damper controller that adjusts the voltage level to MR damper 
thus to track the required damping force. For the system controller a model reference skyhook Sliding Mode 
Controller (SMC) is employed and for damper controller a continuous state algorithm is constructed to determine the 
input voltage thus to gain the required damping force. Along with the proposed modified algebraic model and a 
suitably designed skyhook SMC with continuous -state damper controller, the MR damper is applied to a quarter-car 
suspension model and therefore the performance of the semi-active controller is compared to the active controller and 
with the prevailing passive suspension system. A computer simulation is executed to prove the effectiveness and 
robustness of the semi-active control approach. 

KEYWORDS: Semi- Active Suspension System, Skyhook Sliding Mode System Controller, Continuous-State Damper 
Controller and Quarter-Car Suspension Model 



A suspension, in its more classical and traditional configuration is constituted by three main components: 

• An elastic component (typically a coil spring), that delivers a force proportional and opposite to the suspension 
elongation; this part carries the complete static load. 

• A damping component (typically a hydraulic shock absorber), that delivers a dissipative force proportional and 
opposite to the elongation speed; this part delivers a negligible force in steady-state, however plays a vital role with in the 
dynamic behaviour of the suspension. 

• A set of mechanical components that links the suspended (sprung) body to the unsprung mass. Roughly speaking, 
the suspension could be a mechanical low-pass filter that attenuates the consequences of a disturbance (e.g. an irregular 
road profile) on an output variable. The output variable is usually the body acceleration when comfort is that the main 
objective; the tire deflection when the design goal is road-holding. For passive suspension from Figure- 1 it's clear that 
these two objectives are somehow conflting the tuning and therefore the design of a mechanical suspension tries to search 
out the most effective compromise between these two goals [1]. In this respect, the birth of electronic suspensions for the 
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car mass-market can probably be dated from the early Sixties, when Citroen introduced hydro -pneumatic suspensions [1], 
which modified opened the way to the concept of "on-line" electronic adaptation of the suspension. Heavy power 
consumption, large and unsafe hydraulic systems, and unsure management of the safety issues: the fatal attraction for fully- 
active electronic suspensions lasted solely a few years. They were banned in Fl competitions with in the early nineties, and 
that they never had a big impact on mass-market automotive production [1]. Within the second half of the nineties, a new 
trend emerged: it became clear that the best compromise of cost and performance lay in another technology namely, semi- 
active suspensions. From the paper by Song and Ahmadian [2], it's clear that the semi-active control system presented is 
stable despite of the damping tuning approach. One among the most option of semi-active control systems is that they are 
fail-safe. This suggests that if the control system fails for any reason (including power failure and sensor failure), the 
system acts as a passive system. Therefore a semi-active suspension systems are getting more popular because they provide 
both the reliability of passive systems and also the versatility of active systems without imposing heavy power demands. 
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Figure 1: Filtering Effect of a Passive Suspension: Example of a Road-to-Chassis Frequency 



SEMI-ACTIVE SUSPENSIONS 

Semi-active suspensions are an incredible mixture of appealing features; like negligible power-demand, safety, 
low cost, low weight and significant impact on the vehicle performance. Many types of vehicle are equipped (or are being 
equipped) with semi-active suspension. Such vehicles vary from small vehicles like motorcycles, snowmobiles, etc. to 
massive off-road vehicles, passing through classic cars, and duty- vehicles like trucks, ambulances, fire-trucks, etc. [3-12]. 
It should particularly bear in mind the actual fact that the most semi-active control challenge lies within the dissipative 
constraint of the damper control and not within the spring nonlinearities description [13]. Hence semi-active suspensions 
are implemented by the use of controllable shock absorbers. The technologies out there are based mostly on devices with 
variable orifices (electrohydraulic dampers) or on devices with fluids capable of varying their viscosity as a function of 
electric or magnetic field (electro-rheological and magneto-rheological dampers). Recently, the semi-active suspension 
based mostly on MR damper has attracted a lot of attention [14-17] due to its quick response characteristic to the applied 
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magnetic fields, in-sensitiveness to temperature variations or impurities within the fluid, obtainment of convenient power 
and broad control bandwidth. However, the practical use of MR dampers for control is considerably hindered by its 
inherently hysteretic and highly nonlinear dynamics. Therefore, a dynamic hysteresis model is required to simulate the 
hysteresis phenomenon of MR dampers. To this end, numerous models are proposed within the literature like parametric 
viscoelastic-plastic model based mostly on the Bingham model [18], the Bouc-Wen model [19], non-parametric models 
[20], and lots of. Additionally, theoretical and experimental researches have demonstrated that the performance of a semi- 
active control system is also highly dependent on the selection of control strategy [21]. Therefore, some semi-active control 
schemes are presented and compared in [22] and lots of other approaches, such as skyhook, ground-hook and hybrid 
control [23], Hoc control [24], model-following sliding mode control [25], Neuro-fuzzy control [26] and observer-based 
control [27] are evaluated in terms of their applicability in practice. Thus MR dampers in semi-active vehicle suspension 
applications are determined by two aspects: one is that the accurate modelling of the MR dampers and the other is that the 
choice of an appropriate control strategy. 

In this paper, the control of the stationary response of a quarter car vehicle model to random road excitation is 
considered with semi-active MR dampers. The MR damper is modelled by the proposed modified algebraic model. Except 
for its accuracy, modified algebraic model is additionally more preferable in terms of its low computational expenses 
compared to differential modified Bouc-Wen' s model that is extremely computationally demanding. The control part 
consists of two nuzzled controllers, one the system controller that generates the required damping force, and other the 
damper controller that adjusts the voltage level to MR damper thus to track the required damping force. For the system 
controller a model reference skyhook Sliding Mode Controller (SMC) is employed and for damper controller a continuous 
state algorithm is constructed to determine the input voltage thus to gain the required damping force. Along with the 
proposed modified algebraic model and a suitably designed skyhook SMC with continuous -state damper controller, the 
MR damper is applied to a quarter-car suspension model and therefore the performance of the semi-active controller is 
compared to the active controller and with the prevailing passive suspension system. 

The rest of this paper is organised as follows. Section 3 outlines the description of the modelling of the MR 
damper. Section 4 gives an overview of the quarter vehicle model. Section 5 describes the semi-active control system. The 
results obtained for random road disturbance inputs are presented and discussed in section 6. 
MODELING THE HYSTERETIC BEHAVIOR OF MR DAMPER 

The MR damper employed in the vehicle model with semi-active suspension is an RD-1005 MR damper 
(madeby Lord Corporation Ltd). It's a twin tube MR damper whose conventional, the actual assembly and the components 
are shown in Figure. 2. 
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Figure 2: The Structure of an RD-1005 Damper [28] 
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As a controllable damper, it's subject to the maximal input displacement of 52 mm, the maximal voltage of 12 
V. RD-1005 MR damper is tested by Guo and Hu [28] for sinusoidal excitation with a stroke length of 15mm and a fixed 
frequency of 2 Hz. The test has been performed for five cycles for voltages of 1.0, 2.0, 3.0, 5.0 and 7V. The measured 
force-velocity data for the RD-1005 MR damper is shown in Figure-3. 
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Figure 3: Force Vs. Velocity of RD-1005 MR Damper at 2 Hz Sinusoidal Excitation [28] 

The algebraic model proposed by Guo and Hu [28] is adopted and modified to provide more correct results. The 
model is given by 

F{t) = fo + C b m + -f y tan" 1 {k[x(t) 

% (1) 

-x 0 sgnO(0)] } 

Where F represents the damping force of the MR damper, f 0 the preload of the nitrogen accumulator, C b the 

coefficient of hysteretic damping, f y the yielding force, k the damper coefficient, 0 , v } and v } are the hysteretic 
velocity, the excitation velocity and acceleration of the piston within the damper, respectively. This mathematical model is 
developed primarily based on some physical phenomena. Whereas the primary term is to represent the preload force of the 
pressurized nitrogen gas within the accumulator, the second term is to explain the viscous force of the damper and the last 
one is to replicate the observed hysteretic behaviour, respectively. The mathematical descriptions of the first two terms 
come back from classical mechanics, whereas of the third one is developed primarily based on the definition of a 
trigonometric arctangent function that best resembles the characteristic force-velocity curve of the damper. Further, the 
two terms within the braces of the arctangent function are to account for the lag within the force response to a sinusoidal 

excitation. In this model x{t) = ^(m)M) = a(Dcos(co?),*(0 = -aco 2 sin(cw) 

CO X 

where a is the displacement amplitude and is that angular velocity. In equation C b ,f y , k, and 0 are 

the unknown parameters and to be determined on the basis of experimental data by using least-square curve fitting 
technique. 

In order to validate the algebraic model, Guo and Hu [28] compared the measured damper force and the 
predicted damper force obtained from the algebraic model are shown in Figure- 4 
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Figure 4: Force Vs. Velocity Comparisons between the Algebraic Model Predictions and Experimental Data 

It is observed that there's a general sensible agreement between the estimated and measured values apart from 
higher voltage inputs of 5V and 7V at the highest excitation velocity of 200mm/s (see also Figure-5). The measured force- 
velocity data for the MR damper presented in Figure-3 suggests nonlinear dependence of the force on the applied voltage. 
Starting from this point, the model given in Eq. (1) is modified by multiplying an incremental nonlinear voltage function in 
order to improve the agreement. 
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The nonlinear incremental behaviour of the voltage is characterised by an asymmetric sigmoid function with a bias 
within the lateral axis [29] . The function should also exhibit post-yield limiting behaviour of the damping force attributed 
to the rheological properties of the MR fluid. The nonlinear voltage function k 2 and a 2 are positive constants and I 0 is an 
arbitrary constant representing the bias. The parameters are determined on the basis of experimental data by using least- 
square curve fitting technique. 

For instance, one arrives at the subsequent mathematical model of RD-1005 MR damper for model parameter estimates, 
when a ~ 15mm and f = ^ = 2Hz , respectively. Substituting m = a ® co < m ) mt0 equation (1) yields 
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hysteretic loops of restoring force with respect to the velocity. 
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hysteretic loops of restoring force with respect to the velocity. 
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Figure 5: Hysteretic Loops of Damping Force of RD-1005 MR Damper with Respect to Velocity Obtained from 

the Model given by Eq. (1) at Voltage Inputs of 5 and 7V 
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Figure 6: Hysteretic Loops of Damping Force of RD-1005 MR Damper with Respect to Velocity Obtained from 

the Model given by Eq. (2) at Voltage Inputs of 5 and 7V 

It is obvious from Figure-6 that the proposed modified algebraic model removed the disagreement at the mention 
higher voltage input higher velocity region. This is often presumably because of the impact of the multiplied incremental 
nonlinear voltage function to the algebraic model. It had been inferred that the suggested modified algebraic model might 
overcome the shortcomings of the original algebraic model. Except for its accuracy, modified algebraic model is also more 
preferable in terms of its low computational expenses compared to differential modified Bouc-Wen's model that is highly 
computationally demanding. It's hoped that the current improved model can aid to develop more practical control 
strategies and algorithms for MR dampers. 
SEMI-ACTIVE MR SUSPENSION SYSTEM 

The schematic of the quarter car model with semi-active suspension provided by MR dampers is shown in Figure-7. 
The MR damper is modelled by the modified algebraic model as within the previous section. The equations of motion of 
the semi- active MR suspension system is given by 

m i'Zi =-C 2 (z l -Z 2 )-k 2 (z l -Z 2 )-k l (z l - q) 

+f d -F r +m x g (5) 
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m 2 z 2 = -c 2 (z 2 -z x )-k 2 (z 2 - z x ) 

-f d +K+rn 2 g (6) 
Where m2 is that the sprung mass of the vehicle body, ml is that the unsprung mass, z2 is the absolute 

displacement of the vehicle body, zl is that the unsprung mass displacements, k2 is that the passive suspension stiffness, c2 

is that the passive suspension damping coefficient, kl is that the tyre stiffness, g is that the are the gravity acceleration, and 

Fr is that the friction of suspension, q is that the road(random) excitations and fd is thatthe control force generated by the 

MR damper. The passive components will guarantee a minimal level of performance and safety, whereas the semi-active 

components will be designed to further improve the performance. 



m 2 g k 2 < C 2 



f: 



t, 



mig 



, then the state-equations and output-equations can be listed as follows: 



Figure 7: Semi-Active Suspension System- Quarter Car with MR Damper 

Defining the state variable 
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(8) 



0 
0 

1 

m l 
1 




1 " 








m 2 




q 






F = 


0 


Q = 


g 


U={f d }j = 


Z x ~q 




0 








_z 2 -z { 


m 2 











THE SEMI-ACTIVE CONTROL SYSTEM USING MR DAMPERS 

The semi-active control system based primarily on an MR damper is illustrated as block diagram in Figure-8. The 
plant is that the two degree-of-freedom semi-active suspension system as explained in previous chapter excluding the MR 
damper. The control part consists of two nuzzled controllers. The system controller generates the required damping force 
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whereas the damper controller adjusts the voltage level to track the required force. Within the block diagram, the plant is 
represented by the state space equation (7) and the MR damper has been modelled as given in section 2. The damping force 
of MR damper island the required damping force generated by the system controller is/ c . 
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System 
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Damper 
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Figure 8: Block Diagram of Semi- Active Control System with MR Damper 

System Controller: Skyhook Sliding Mode Controller 

This section provides an summary of the system controller, that was derived in reference [30]. In designing a 
Skyhook Sliding Mode Controller (SMC), the target is to consider the subsequent nth order non-linear system because the 
controlled plant, and thus defined by the general state-space within the equation 



(9) 



x = f(x,u,t) 

Where xER n is that the state vector, n is that the order of the nonlinear system, and u G R m is that the control input, m 
is that the number of inputs. 

A time varying surface s(t) is outlined within the state space R (n) . S (e, t) is that the sliding surface of the hyper- 
plane, which is given in equation (10) and shown in Figure-9. 



s(e,t) = \^— + X 

In the 2-DOF active suspension system, we let n = 2, given that, because it is a second-order system. 
s = e + Xe 



(10) 
(11) 



sliding surface 
Slope is X 




Figure 9: Sliding Surface Design 

Where velocity and acceleration of the vehicle body are selected as error (e) and change in error ( ) feedback 
signals for the 2-DOF SA suspension system control and X is a positive constant that defines the slope of the sliding 
surface. 

From the above equation, the second-order tracking problem is currently being replaced by a first-order 
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stabilization problem within which the scalar s is kept at zero by means of governing condition [31]. 

V(s) = U (12) 

This is obtained from the utiliation of Lyapunov stability theorem and it states that the origin is a globally 
asymptotically stable equilibrium point for the control system. Equation shown on the top of is positive definite and its 
time derivative is given in inequality, 



(13) 



V(s) = ss<0 

To satisfy the negative definite condition, the system should satisfy the above shown inequality. If the inequality 
equation is not satisfied means, the state of system deviates from the sliding surface. It is desirable that whenever the 
system isn't on the sliding surface, at a given instant of time it's driven on the sliding surface and kept there. The variable s 
will be driven to zero by means that actual control input u from the controller. 

The controller design procedure consists of two steps. First, a feedback control law u is chosen to verify sliding 
condition. However, so as to account for the presence of modelling imprecision and of disturbances, the control law needs 
to be discontinuous across the s(t). If the implementation of the associated control switching is imperfect, this may lead to 
chattering. Thus, in a second step, the discontinuous control law u is suitably smoothed to attain an optimal trade-off 
between control bandwidth and tracking precision. The skyhook control will reduce the resonant peak of the sprung mass 
quite considerably and therefore can achieve a good ride quality. By borrowing this concept to reduce the sliding chattering 
phenomenon, as in Figure- 10 shown, a soft switching control law is introduced to the major sliding surface switching 
activity of control law that is to attain a sensible switch quality for the SkyhookSMC. 
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Figure 10: Sliding Mode Surface Design with Skyhook Control Law 

This can be done by stipulating the following condition. 



= fc=\ 



-c n tanh 



ss >0 



(14) 



0 ss < 0 

The variable of s is outlined in Equation (10), that contains the system information. It may be taken as a part of the 
Skyhook SMC control law in Equation (13), where c 0 is an assumed positive damping ratio for the switching control law. 
The Skyhook SMC has to be chosen in such a way that the existence and also the reachability of the sliding-mode are both 
guaranteed. Noting that S is an assumed positive constant that defines the thickness of the sliding mode boundary layer [32]. 
The control system parameters employed in this research are identical as those used in [32] and are summarized in Table 1. 



Table 1: System Controller Parameters 
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Parameter 


Symbol 


Value 


SkyhookSMC damping 
coefficient 


c 0 


-5000 


Thickness of the sliding 
mode boundary layer 


8 


28.1569 


Sliding surface slope 


X 


10.6341 



Damper Controller: Continuous-State Control 

Force-tracking control of the MR damper model is performed employing a simple continuous-state control. The 
continuous-state control of the MR damper was employed by [33, 34 and 35], within which a simple feedback control 
strategy might linearize the response of the MR damper. To induce the MR damper to give roughly the corresponding 
required control force, the command signal is chosen and described by the SIMULINK block diagram as in Figure- 1 1. 




Figure 11: Force Tracking Control of MR Damper 

The damping force of MR damper is fed back with a gain B and compared to the required force. The resultant 
error is scaled by a gain K. To make sure the damper cannot generate energy to the system, the controller function is 
enabled solely when the direction of damping force and also the error have the same direction. Therefore it's essential to 
have a sign correction. If they have different sign, the input voltage has to set to be zero. The acceptable control signal is 
then varied between the maximum and minimum voltages. The continual state algorithm for selecting the input voltage 
may be stated as 



if 

G(f c - Bf d ) sgn(/ d )>V.^v = 
elseif 

G(f c - Bf d ) sgntf, ) < -+ v = 
else 

v = G(f c -Bf d ) Sga (f d ) 
Where, 

^ max is the voltage to the current driver associated with saturation of the MR damper 

^ min is the minimum voltage to the damper (i.e. 0V) 
f 

J c is the required control force determined by the system controller 



(15) 
(16) 
(17) 



is the damping force of the MR damper 



The values of G and H are decided by the trial and error method. In this paper, the values of K , H and V max were set 
to be 0.0038 V/N , 1, and 2 V respectively, as in [36]. 



RESULTS AND DISCUSSIONS 
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To evaluate the performance of the semi-active suspension controller with MR damper, three kinds of suspensions, 
particularly passive, active and semi-active suspension, are studied in this work. The quarter car vehicle model with the 
subsequent parameters is considered as an example. m2=240kg, m 1 =3 Okg, k2= 1 6kN/m, c2=1400N-s/m, kl = 16()kN/m, g=9. 
81m/s2 and Fr=300N. In the simulation of the control on two degrees of freedom suspension system, it's excited by a 
random road disturbance loading that is described by the road profile with the parameters of reference space frequency 
n0=0. lm-1, road roughness coefficient P (nO) =256*10-6m3/cycle and vehicle speed v0=72km/h. The velocity and 
acceleration of the vehicle body are selected as error (e) and change in error (EC) feedback signals for the 2-DOF active 
suspension system control. In so doing, it's expected to lead to reduced values of sprung mass acceleration, suspension 
operating space and dynamic tire load, and hence, improved suspension performance. 

Passive suspension means that the control force is set to zero for all time, the passive damping is equal to 1400 
Ns/m. Active suspension means control input fc is totally realised by Equation (13). Semi-active suspension means the 
control input fd is realised by the MR damper with the control structure in Figure 8 and 1 1 . The input voltage of the MR 
damper is restricted inside 0 ~ 2V. The numerical results are obtained using a MATLAB/SIMULINK. 

There are three performance indexes for the vehicle suspension system, which includes body acceleration, 
suspension deformation and tire deformation. In this context, the results for the three indexes are applied to evaluate the 
performance for the ride comfort of the 2-DOF semi-active vehicle suspension system. 

Figure - 12, 13 and 14 illustrates clearly how semi-active controller primarily based on the proposed modified 
algebraic model will effectively absorb the vehicle vibration as compared to the active and also the passive system. Figure. 
12 shows the sprung mass acceleration of all active, semi- active and passive suspension system for comparison purposes. 
The body acceleration within the semi-active suspension system is reduced considerably, that guarantee better ride comfort. 



1/4 vehicle suspension system body acceleration 




□ 1 23456789 10 
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Figure 12: Sprung Mass Acceleration Response 

Figure. 12 and 13 shows the suspension deflection and tire load of all active, semi- active and passive suspension 
system for comparison purposes. The suspension deflection and tyre load is also smaller using the semi-active controller 
based on the proposed modified algebraic model, that guarantee better road holding. 
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1/4 vehicle suspension system suspension distorsion 
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Figure 13: Suspension Stroke Response 
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Figure 14: Tyre Load Response 



Thus in overall from Figure. (12- 14), it may be seen that both active and semi-active suspension systems are able to 
comparatively lower magnitude for vehicle body acceleration, suspension deflection and tire load within the time domain 
in comparison with the passive suspension system. Using the presented control structure in Figure-8, the semi-active 
suspension system in conjunction with the MR damper can achieve a control performance that's almost like that of the 
active suspension system apart from a little deterioration due to the passivity and the limitation constraints. It demonstrates 
the effectiveness of the semi-active controller with MR damper for vibration suppression of the suspension system. 

Figure. 15 shows the psd(power spectral density) of sprung mass acceleration for various suspension systems. The 
vehicle model considered for the analysis is two degrees of freedom vehicle model, thus two predominant frequencies are 
observed within the Figure-15. It may be seen that, the sprung mass acceleration characteristics as compared to the passive 
suspension, is improved by the active and semi-active suspension system at both the predominant frequencies. 
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Figure 15: Sprung Mass Acceleration Psd 

All these results indicate that the semi-active controller presented in this paper dissipate the energy due to road 
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excitation very well and improve both the ride comfort and vehicle stability. 
CONCLUSIONS 

The paper presents a new modified algebraic model for semi-active MR damper suspensions in controlling the 
stationary response of a quarter car vehicle model traversing a rough road with constant velocity. The modified algebraic 
model will accurately characterize the dynamics of the MR damper. A static output skyhook sliding mode controller is 
designed directly using the measurable sprung mass velocity and acceleration signals that generate the required damping 
force. Next a continual state algorithm is constructed to determine the input voltage so as to gain the required damping 
force. Along with the proposed modified algebraic model and a suitably designed skyhook sliding mode controller with 
continuous-state damper controller, the MR damper is applied to a quarter-car suspension model. The performance 
characteristics and the robustness of the semi-active suspension system are evaluated by two nuzzled controllers, and then 
compared with the active and passive suspension system. The results show that the performance of semi-active controller 
primarily based on the proposed modified algebraic model is better than the performance of the passive suspension and can 
achieve compatible performance as that of active suspension controller apart from a little deterioration. 
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